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1 Introduction 

Doubtlessly the most central notion in superconductivity, for both low and 
high transition temperatures Tc, is that of Cooper pairs (CPs) [1] that form 
among the underlying electron (or hole) charge carriers of the many-electron 
system. And yet, it is perhaps the least understood concept. Shortly after 
the publication of the BCS theory [2] of superconductivity, charged pairs 
observed in magnetic flux quantization experiments with 3D conventional 
[3] [4], and much later with quasi-2D cuprate [5] superconductors, suggested 
CPs as an indispensable ingredient, regardless of the inability of BCS theory 
per se to describe high-Tc superconductivity, as now seems an almost uni- 
versal consensus. Cooper pairing is no less central, albeit with a different 
interfermion interaction, to neutral-atom superfluidity as in liquid ^He [6], 
and presumably also to ultracold trapped alkali Fermi gases such as ^Li [7] 
and '^'^K [8] where pairing is expected to occur as well. More recently, stud- 
ies of quantum degenerate Fermi gases consisting of neutral "^^K atoms and 
their so-called Feshbach "resonance superfluidity" have appeared [9]- [13], but 
assuming quadratically-dispersive CPs and ignoring hole pairs — two funda- 
mental shortcomings as we now illustrate. 

In Sec. 2 we recall how the Bethe-Salpeter (BS) many-body equation (in 
the ladder approximation), treating both 2p and 2h pairs on an equal footing, 
implies that the ordinary CP problem [based on an ideal Fermi gas (IFG) 
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ground state (the usual "Fermi sea")] does not possess stable energy solu- 
tions; in Sec. 3 we sketch how CPs based not on the IFG-sea but on the BCS 
ground state survive, along with the usual trivial sound mode, as nontrivial 
"generalized" or "moving" CPs, linear in total or center-of-mass-momentum 
(CMM) in leading order, that are positive energy resonances with an imag- 
inary energy term implying finite-lifetime effects. The nontrivial "moving 
CP" solution, though often confused [14] with it, is pliysically distinct from 
the trivial sound mode solution sometimes called the Anderson- Bogoliubov- 
Higgs (ABH) [15], ([16] p. 44), [17] [18] collective excitation mode. The ABH 
mode is also linear in CMM in leading order, and reduces to the IFG ordinary 
sound mode in zero coupling. All this occurs in both the 3D study outlined 
in Ref. [19] as well as in 2D [20]. Sec. 4 offers conclusions. 

Provided CPs are bosons, as they indeed [21] turn out to be, our results 
will in general be crucial for Bose-Einstein condensation (BEC) scenarios em- 
ploying BE models of superconductivity, not only in exactly 2D as with the 
Berezinskii-Kosterlitz-Thouless [22] [23] transition, but also down to (1 + e)D 
which characterize the quasi-lD organo-metallic (Bechgaard salt) supercon- 
ductors [24]- [26]. This is contrary to well-entrenched perceptions (see, e.g., 
Ref. [27]) that BEC is impossible in 2D. 



2 Ordinary Cooper pairing 

The original "ordinary" CP problem [1] (with the BCS model two-electron 
interaction) is defined in an A^-electron system for two electrons above the 
Eermi energy Ep but within a thin shell of energy hwo- Only there do they 
suffer a constant attraction —V < 0, giving rise to the familiar result for the 
negative energy of the pair, relative to the energy 2Ef of the interactionless 
pair, 

^o = -4^ — -2^^e-/\ (1) 

Here A = g{EF)V is a dimensionless coupling constant with glEp) the density 
of fermionic states (for each spin) evaluated at Ep. The equality in (1) is 
exact in 2D for all coupling — as well as in ID or 3D provided only that 
<^ Ep so that g{e) ~ g[Ep), a constant that can be taken outside an energy 
integral. 
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However, the original CP problem neglects the effect of two-hole (2h) CPs 
treated simultaneously on an equal footing with two-electron, or two-particle 
(2p), CPs — as Green's functions [28] can naturally guarantee. On the other 
hand, the BCS condensate consists of equal numbers of 2p and 2h Cooper 
correlafAons. This was already evident, though scarcely emphasized, from 
the perfect symmetry about fi, the electron chemical potential, of the well- 
known Bogoliubov [29] v'^{e) and M^(e) coefficients [see just below (10) later 
on] , where e is the electron energy. Indeed, our prime motivation comes from 
the fact established recently [30] that the BCS condensate is a BEC conden- 
sate for equal numbers of 2p and 2h pairs, in the limit of weak coupling. 
Additional empirical motivation comes from the unique but unexplained role 
played by hole charge carriers in the normal state of superconductors in gen- 
eral [31]. Even further motivation stems from the ability of the "complete (in 
that both 2h- and 2p-CPs are allowed in varying proportions) BE model" of 
Refs. [30], [32]- [33] to "unify" both BCS and BEC theories as special cases, 
and to predict substantially higher Tc's than BCS theory without abandoning 
electron-phonon dynamics. Compelling evidence for a significant presence of 
this dynamics in high-Tc cuprate superconductors from angle-resolved pho- 
toemission spectroscopy data has recently been reported [34]. 

In dealing with the many-electron system we assume the BCS model 
interaction in the form with double Eourier transform 

i/(|ki - k;|) = -(kl/kik[)V if kF-kD< ki, k[ <kF + ko, (2) 

and = otherwise. As before V > 0, hkp = mvp the Eermi momentum, m 
the effective electron mass, vp the Eermi velocity, and ko = ojd/vf with ujd 
the Debye frequency. The usual physical constraint fiwo ^ Ep = h'^kp/2m 
then implies that k^/kp = liwr)/2Ep <C 1. 

The bound-state BS wavefunction equation [19] in the ladder approxima- 
tion with both particles and holes for the original lEG-based CP problem 
using an interaction such as (2) is 
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where L'^ is the "volume" of the rf-dimensional system; K = ki + k2 is the 
CMM and k = |(k_^ — k2) the relative wavevcctors of the 2e bound state 
whose wavefunction is ^(k,£'); = Ei + E2 is the energy of this bound 
state while E = Ei — E2, and Gq (K/2 + k, £^/2 + is the bare one-fermion 
Green's function given by 



i [ -El + eki - Ep-ie -Ei + ekj - Ep + ie 



(4) 



where Cki = TT?k\/2m and 9{x) = 1 for a; > and = for a; < 0, so that 
the first term refers to electrons and the second to holes. Figure 1 shows 
all Feynman diagrams for the 2p, 2h and ph wavefunction ip^, ip- and ipQ, 
respectively, that emerge in the general (BCS-ground-state-based) problem 
to be discussed later. For the present IFG-based case, diagrams in shaded 
rectangles do not contribute. Since the energy dependence of ^(k, E) in (3) 
is only through the Green's functions, the ensuing energy integrals can be 
evaluated directly in the complex i5'-plane and yield, for interaction (2), 

(2a - £:o) = F J] Vk' -VY, 'Vk' (5) 

k' k' 

where ■^k is the resulting wavefunction after the energy integration. Here 
a = fv^k"^ /2m — Ep while £q is the (unknown) cigcnvahic energy. The single 
prime over the first (2p-CP) summation term denotes the restriction < 
< hwo (i.e., above the IFG "sea") while the double prime in the last 
(2h-CP) term means —hwo < Cfe' < (i.e., helow the IFG sea). Without this 
latter term we have Cooper's Schrodinger-hke equation [1] for 2p-CPs whose 
implicit wavefuncion solution is clearly ■^k = (2^ — ^^o)"''^^ X^k' ^k'- Since 
the summation term is constant, performing that summation on both sides 
allows canceling the ■?/'k-dependent terms, leaving the eigenvalue equation 
Sk(2Cfe~^o)~^ = This is one equation in one unknown £q] transforming 

the sum to an integral over energies immediately gives (1). This corresponds 
to the usual negative-energy, infinite-lifetime stationary-state bound pair. 
For 7^ ^ the CP eigenvalue equation is just (2^ — £'i^ + /i^X^/4m)~^ = 
Since a CP state of energy £k is characterized only by a definite K but 
not definite k, in contrast to a "BCS pair" defined [Ref. [2], Eqs. (2.11) 
to (2.13)] with fixed K and k (or equivalently definite ki and k2). This 
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renders CPs legitimate bosons [21]. Without the first summation term in (5) 
the same result in Sq (1) for 2p-CPs follows for 2h-CPs (apart from a sign 
change) . 

In contrast with Cooper's equation neglecting hole-pairs, the complete CP 
equation (5) cannot be derived from an ordinary (non-BS) Schrodinger-like 
equation in spite of its simple appearance. A more general technique such 
as the BS equation that includes both particles (in this case electrons) and 
holes is needed. To solve it for the unknown energy £o, let the rhs of (5) be 
defined as A — B, with A relating to the 2p-pair term and B to the 2h-pair 
term. Then the unknown i/j^ becomes 

V'k = (^ - B)/{2^k - So) or equivalents ^(0 = {A - B)/{2i - £o) 

(6) 



whence 



A = X = oi^-B)X dz/z ={A- B)x, 

Jo ^ J-Eo 

B = xf = 1{A -B)X [ ° dz/z = iA- B)y. (7) 

The integrals are readily evaluated giving x = |Aln(l — 2hujD/So) and y = 
— |Aln(l + 2hj£)/So). As A and B still contain the unknown let us 

eliminate them. Note that equations (7) are equivalent to two equations in 
two unknowns A and B, or 

{l-x)A + xB ^0 and - yA + {1 + y)B ^ 0. 

These readily lead to the single equation 1 — x + y = 0, which on inserting 
the definitions for x and y becomes 

l^-Xln[l-{2hu}D/Sof] which gives So = ±%2hh)BN e^l^ - 1- 
2 (8) 

As the CP energy is pure-imaginary, there is an obvious instability of the CP 
problem when both particle- and hole-pairs are included. This transcendant 
result dates back to the late 50's and early 60's and was reported in Refs. [16] 
p. 44 and [35] Sec. 33, where, however, the well-known pure 2p and 2h special 
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cases just stressed was not discussed. Clearly then, the original CP picture 
is meaningless if particle- and hole-pairs are treated on an equal footing, as 
consistency demands. Curiously, this result has been largely ignored in the 
entire literature since then. 

A dramatic analogue of the nontriviality of such consistency is found in 
the very high temperature treatment of relativistic BEC [36], where pair pro- 
duction becomes possible and creates antibosons in addition to more bosons. 
Here, BEC must take into account N antibosons of charge, say, —q along 
with the N bosons of charge q. In units such that h = c = ks = 1 the 
boson energy is €k = {K'^ + m^)^/^. Charge conservation requires that not 
only N = iVo(T) + '^j^^qI'SW P{^k — I^b) — be constant but rather 
A^— A^, where A^ is the same as A^ but with instead of —/is- If 

p = g(A^ — N)/L^ = qn is the net conserved charge density, it is shown 
in Ref. [36] that = {3\n\/mB)^ and that the condensate fraction rio/n = 
[1 — (T/Tc)^]. This is qualitatively different from the better- known results as- 
suming only A^ constant, which are the mass- independent Tc = [n'^n/(^{3)Y^^ 
and no/n = [1 — (T/Tc)^]. This example exhibits the strikingly dramatic ef- 
fect of including or not antiparticles (analogous to holes in the nonrelativistic 
case) . 

3 Generalized Cooper pairing 

However, a BS treatment not about the IFG sea but about the BCS ground 
state (which we refer to as "generalized" Cooper pairing) vindicates the CP 
concept, and adds something new. This is equivalent to starting not from the 
IFG unperturbed Hamiltonian but from the BCS one. Thus, (4) is replaced 

by 

where = ^/^]^~+~^ with A the fermionic gap, vl = |(1 — Ck/E^) and 
= 1—vl are the Bogoliubov functions [29]. As A — >■ these three quantities 
become |,^^.|, 9{ki — kp) and 9{kp — ki), respectively, making (9) become (4) as 
expected. Substituting Go(ki, Ei) by Go(ki, Ei) corresponds to rewriting the 
total Hamiltonian so that the pure-kinetic-energy unperturbed Hamiltonian 
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is replaced by the BCS one. The remaining terms are then assumed suitable 
to a perturbation treatment. Experimental support for this can be found 
precisely in Refs. [3]-[5], and its physical justification lies in recovering both 
the expected ABH sound mode (which contains the BCS T = gap equation) 
and the finite-lifetime effects of moving CPs. In either 3D [19] or 2D [20] 
the BCS-based BS equation yields two distinct solutions: a) the usual trivial 
ABH sound solution and b) a highly nontrivial "moving CP" solution. In 
either case the BS formalism consists of a set of three coupled equations, 
one for each (2p, 2h and ph) channel wavefunction for any spin-independent 
interaction such as (2). However, the ph channel decouples, leaving only 
two coupled wavefunction equations for the ABH solution in 2D which we 
consider first. We focus here on 2D because of its interest [37] for quasi-2D 
cuprate superconductors. 

The equations involved are too lengthy even in 2D and will be derived in 
detail elsewhere, but for the trivial or ABH sound solution they boil down 
to the single expression 

— XhVF / dk d(p{uK/2+]<:UK/2-k + t'K/2+k l^K/2-k} X 



X 



' kp—ko 

■f^K/2+k'J^K/2-k _^ MK/2+k«K/2-k 



+ -E'K/2+k + -E'K/2-k —^K + -E'K/2+k + -E'K/2-k, 



= 1 

(10) 



where (/? is the angle between K and k. Here ki) is defined as below (2) and 
as before A = Vg{EF) with g^Ep) = m/27ih'^ the constant 2D electronic 
DOS and V is defined in (2). The ABH collective excitation mode energy 
Sk must then be extracted from this equation. For K = it is just Sq = 
(Ref. [16] p. 39) and (10) rewritten as an integral over ^ = h'^k^/2m — Ep 
reduces to the familiar BCS T = gap equation /q^" rf,^/ a/.^^ + = 1/A 
for interaction (2) which gives A = tiwi)/ sinh(l/A). Returning to the ABH 
energy £k equation (10) and Taylor-expanding £k about K — and then 
taking A small, leaves 

= + 0{K^) + 0(A) (11) 

where o(A) refers to interfermion interaction terms that vanish as A — 0. 
Note that the leading term is just the ordinary sound mode in an IPG whose 
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sound speed c = vpl\fd, in d dimensions. This result also follows (trivially) on 
solving for c in the familiar thermodynamic relation dPj dn = nic^ involving 
the zero-temperature IFG pressure P — 'n?[d{E / N) / dn\ — 2nEF/{d + 2) = 
2Cd'n?/'^ '^^ I (d-|-2) where the constant will drop out. Here the IFG ground- 
state energy E = dNEp/{d + 2) was used along with Ep = h'^kj^/2m = 
C^n^/"^ while n = N/L'^ = kp/ d2'^~'^n'^/'^V{d/2) is the fermion- number density. 
The derivative dP/dn finally gives c = hkp/m\fd = vp/y/d which in 2D is 
just the coefficient of the first term of (11). 

The nontrivial or moving CP solution of the BCS-ground-state-based BS 
treatment, which is entirely new, leads to the pair energy Sk which in 2D is 
contained in the equation 

j-kp+ko j-2-K 

— XhVp / dk I rf97MK/2+kt'K/2-k X 

^"^ Jkp-kn Jo 

r -I -E'K/2+k + -E'K/2-k 

X |MK/2-k^^K/2+k - «K/2+k^K/2-k) p2 i / f ^2 = 

-^K + (-C'K/2+k + -C/K/2-kj ^12) 

In addition to the pp and hh wavefunctions (depicted diagrammatically in 
Ref. [19] Fig. 2), diagrams associated with the ph channel give zero contri- 
bution at T = 0. A third equation for the ph wavefunction describes the ph 
bound state but turns out to depend only on the pp and hh wavefunctions. 
Taylor- expanding in (12) in powers of K around K = 0, and introduc- 
ing a possible damping factor by adding an imaginary term —iVx in the 
denominator, yields to order K'^ for small A 



±£k ^ 2A + —hvpK + -^e'/^K^ - ihvpK 
2ti 9 kp) 



- + -^e'/'K 
n 12kD 



+ 0{K^) 
(13) 



where the upper and lower sign refers to 2p- and 2h-CPs, respectively. A 

linear dispersion in leading order again appears, but now associated with the 
bosonic moving CP. Figure 2a graphs the exact moving CP real energy (full 
curves) extracted from (12), along with its leading linear-dispersion term 
(short-dashed) and this plus the next (quadratic) term (long-dashed) from 
(13). The interaction parameter values used for (2) were hoJo/Ep — 0.05 (a 
typical value for cuprates) and the two values A = | (lower set of curves) 
and I (upper set), so that E^/Ep = 2/\/Ep = 2^/5/i?Fsinh(l/A).This has 
values ~ 0.004 (for X = ^) and 0.028 (for A = |), (marked as dots in the 
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figure). Remarkably enough, the hnear approximation (short-dashed hues in 
figure) is better over a wider range of Kjkp values for weaker coupling in 
spite of a larger and larger partial contribution from the quadratic term in 
(13). This peculiarity also emerged from the ordinary CP treatment with 
the BCS model interaction [38] in both 2D and 3D, and with a (regularized) 
delta interaction in 2D [39] and in 3D [40], and might suggest the expansion 
in powers of K to be an asymptotic series that should be truncated after the 
linear term. For reference we also plot the linear term hvpK/ -\/2 of the sound 
solution (11) (thick long-dashed line starting at origin). The posifAve-cnergy 
2p-CP resonance of width has a lifetime from (13) of = fi/2Vx = 
h/2 [{\/'K)hvFK + {hvF/12kD)e^/^K'^] diverging only at A' = 0, and falling 
to zero as K increases; see Fig. 2b. Thus, "faster" moving CPs are shorter- 
lived and eventually break up, while "non-moving" ones are (infinite-lifetime) 
stationary states. The real linear term {\/2T{)hvpK in (13) contrasts sharply 
with the coupling-independentlcsiding-teYm {2/7r)hvpK that follows [38] from 
the original CP problem neglecting holes, which if graphed in Fig. 2a would 
almost coincide with the ABH term hvFK/\/2 and have a slope about 90% 
smaller. Fig. 2c depicts analogies of a 3D potential problem for the original 
CP energy (top), as well as for X > (middle) and K — (bottom) BCS- 
based BS CPs. 

As in Cooper's [1] original equation, our BCS-based BS moving CPs are 
characterized by a definite K and not also by definite k as the pairs dis- 
cussed by BCS [2]. Hence, the objection does not apply that CPs are not 
bosons because BCS pairs with definite K and k (or equivalently definite 
ki and k2) have creation/annihilation operators that do not obey the usual 
Bosc commutation relations [Ref. [2], Eqs. (2.11) to (2.13)]. In fact, (12) 
shows that a given "generalized" CP state labeled by either K or £k can 
accommodate (in the thermodynamic limit) an indefinitely many possible 
BCS pairs with different k's; see Ref. [21]. A recent electronic analog [41] of 
the Hanbury Brown- Twiss photon-effect experiment suggests electron pairs 
to be definitely bosons. 

4 Conclusions 

Hole pairs treated on a par with electron pairs play a vital role in determining 
the precise nature of nontrivial CPs even at zero temperature — only when 
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based not on the usual IFG "sea" but on the BCS ground state. Their 
treatment with a Bethe-Salpeter equation gives purely-imaginary-energy CPs 
when based on the IFG, and when based on the BCS ground state gives 
positive-energy, resonant-state CPs with a finite lifetime for nonzero CMM. 
This is instead of the more familiar negative-energy stationary states of the 
original IFG-based CP problem that neglects holes, as sketched just below 
(5). The BS "moving-CP" dispersion relation (13), on the other hand, is 
gapped by twice the BCS energy gap, followed by a linear leading term 
in the CMM expansion about K = 0. This linearity is distinct from the 
better-known but trivial one (11) associated with the sound or ABH collective 
excitation mode whose energy vanishes at i^' = 0. 

Thus, instead of the quadratic h^K"^ /2(2m) assumed for CPs in Refs. [9]- 
[13], [30], [32], [42]- [44], among many others, BF models based on the correct 
CP linearity for the boson component can give BEG for all ci > 1, including 
exactly 2D, and thus in principle address not only quasi-2D cuprate but also 
quasi- ID organo-metallic superconductors. 
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Figure captions. 

1. Wavefunction Feynman diagrams for 2p 2h and ph (iPq) 
bound states arising from the BS equations. Shaded rectangles designate 
diagrams that do not contribute in the IFG-based case. 

2. a) Exact "moving" 2p-CP (real) energy £k (in units of Ep) in 2D from 
(12) (full curves), compared with its linear leading term (short-dashed 
lines) and its linear plus quadratic expansion (long-dashcd curves) both 
from (13), vs. CMM wavenumber K (in units of kp), for interaction (2) 
parameters A = | (lower set of curves) and | (upper set of curves), and 
hwo/ Ep = 0.05. For reference, leading linear term (11) of trivial ABH 
sound mode is also plotted (lower thick dashed line), b) 2p-CP hfetime 
as defined in text, c) Analogy of BCS-based BS 2p-CPs with various 
states in a 3D potential problem, as discussed in text. 
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Figure 1. Wavefunction Feynman diagrams for 2p 2h {ip-) and pli {il^o) bound 

states arising from the BS equations. Shaded rectangles designate diagrams that do 
not contribute in the IFG-based case. 



